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ABSTRACT 

By analysing K2 short-cadence observations we detect starspots on WASP-85A, the host star of the hot 
Jupiter WASP-85Ab. The detection of recurring starspot occupation events indicates that the planet’s 
orbit is aligned with the star’s rotational axis (A < 14°) and suggests a stellar rotational period of 
15.1 =b 0.6 d. The K2 lightcurve reveals a rotational modulation with a period of 13.6 zb 1.6 d, consistent 
with the period determined from starspots. There are no significant transit-timing variations and thus 
no evidence of any additional planet in the system. Given the pronounced rotational modulation we 
are only able to place an upper limit of 100 parts per million for any phase-curve modulations and 
the secondary eclipse. 

Keywords: planets and satellites: individual (WASP-85Ab) - stars: individual (WASP-85A) - 
starspots 


1. INTRODUCTION 


Since the failure of a second reaction wheel, the Kepler 
telescope provides the communi ty wit h K2 observa tions 
of fields along the ecliptic plane (Howell et al. 2014), and 
is therefore observing some previously known WASP ex¬ 
oplanets. One result has been the discovery of two addi¬ 


tional transiting planets in the WASP-47 system (Becker 
eFaLl|2015 ). 


The tra nsit ing hot Jupiter WASP-85Ab was discovered 
by Brown et al. (2014) and was observed by K2 in Cam¬ 
paign 1. It has a 2.66-d orbit around its G5, V = 11.2 
host star. WASP-85 A forms a close visual binary, with an 
angular separation of 1.5arcsec, with a cooler and dim¬ 
mer K0, V = 11.9 WASP-85B. The WASP photometry 
showed that one of the two binary components exhibits 
a 14.6 zb 1.5 d r otationa l modulation, which indicates stel¬ 
lar activity ([Brown et al.||2014). 

When a starspot is occulted by a transiting planet the 
lightc urve e xhibits a “bump” of temporary brightening 
(Silva 2003). The detection of starspots can provide a 
measurement or a constraint on the stellar obliquity, i.e. 
the angle between the stellar rotation and the planet’s 
orbital axis. For example, the recurrence of the same 
starspot in consecutive tra nsits has re vealed that Kepler- 
17 has a low obliquity (Desert et al. 2011). The ability 
to detect the recurring starspot occultation in a partic¬ 
ular system depends on the ratio between the stellar ro¬ 
tational and planet orbital periods, and the lifetime of 
starspots. Alternatively, if starspots are occulted at cer¬ 
tain preferential phases of the transit, this may constrain 
the obliquity if the transit chord crosses the active stellar 
latitudes at t he o b served active phases, a s in the case of 
HAT-P-11 (Sanchis-Ojeda & Winn 201l|). 

The majority of hot Jupiters arouncTstars cooler than 
around 6250 K are found in aligned orbits, suggestive of 
a disc migration mechanism. Hotter stars, on the other 
hand, host hot Ju piters that are in many cases found in 
misaligned orbits (Albrecht et al. 2012), which cannot 


be accounted for solely by planet disc migration mecha¬ 


nism. Instead, it is believed that planet-planet scatter¬ 
ing, Kozai mechanism and tidal dissipation play an ad- 
ditio nal role i n shaping misaligned orbits of hot Jupiters 
(e.g. Triaud et al.||2010). 

In this paper, we present a detection of starspots in the 
K2 lightcurve of WASP-85, which we use to constrain the 
stellar obliquity. We also refine the system parameters 
and the stellar rotational period, and search for transit¬ 
timing variations, as could be caused by an additional 
body in the system. 

2. THE K2 OBSERVATIONS 

The K2 Campaign 1 observations cover an 82-d time- 
span between 2014 May 30 and 2014 August 20 and 
provide a total of 120 660 short-cadence images. We 
retrieved the 1-min short-cadence target pixel file for 
WASP-85 via the Minkulski Archive for Space Telescopes 
(MAST). 

Given the degraded pointing accuracy of K2, the pho¬ 
tometric precision is reduced because the position of an 
observed object drifts within the photometric aperture, 
and because of the n on-unifo rm pixel sensitivity of the 


on-board detectors. Vanderburg & Johnson (2014) de¬ 


veloped a self-flat-fielding (SFF) procedure with which 
they correct the long-cadence photometry to within a 
factor of two of the original Kepler precision by correlat¬ 
ing the measured flux and arclength of the drift. Our K2 
data-reduction procedure for short-caden ce observations 
is similar to th e procedure presented by |Vanderburg fc] 
Johnson] (|2014|). 

We first applied a standard data reduction procedure 


using PyRAF to ols for Kepler (PyKE) version 2.6.2 (Still 
& Barclay 2012). Besides removing all the frames that 
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were originally quality-flagged by K2, we manually ad¬ 
justed quality flags and removed an additional 100 frames 
between 2014 June 1 18:04 UT and 19:41 UT, since they 
were obtained while the spacecraft was still in coarse 
pointing after the attitude tweak. Next, we used a 
large 100-pixels aperture mask to accommodate the di¬ 
luted photometric signal from both binary components, 
WASP-85 A and WASP-85B, and to reduce the photo- 
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Figure 1. SFF calibration procedure. All the panels correspond to a time interval of 5d. la: Position of the center of the object given 
as CCD row vs. CCD column. Spatial outliers are marked with red and flagged as thruster events, lb: Measured normalized flux vs. 
arclength of the spacecraft drift. The correlation is fitted with a Gaussian convolution (marked with red), lc: Lightcurve before the SFF 
correction with the quality-flagged data points removed. Id: Lightcurve after the SFF correction with the quality-flagged data points 
removed, including the thruster events which were identified during the SFF procedure. The positive flux excursions in the corrected 
lightcurve are the remaining cosmic rays that have avoided being identified and quality-flagged as cosmic rays in the downloaded K2 target 
pixel file. 


metric losses due to the relative drift between the object 
and the extraction aperture. Photometry of the individ¬ 
ual binary components could not be performed because 
the angular separation (1.5arcsec) is smaller than the 
Kepler point spread function (full width half maximum 
of 6 arcsec). 

After the data extraction and removal of low-frequency 
modulations, we performed our SFF procedure. The pro¬ 
cedure is based on the PyKE tool kepsff with two main 
modifications to adjust and improve the efficiency for 
the short-cadence data. First, the 1-min short-cadence 
frames obtained during the thruster firing events are seen 
to have object positions that are outliers from the mean 
position drift on the detector. Data points with posi¬ 
tional deviations of more than 6a were flagged as thruster 
firing events and deleted (see Figure la). This level was 
established by trial and error, as rejecting the thruster 
firing events while removing a minimal number of good 
data points. A total of 490 data points were quality- 
flagged as thruster firing events, which represents 0.4 per 
cent of the entire lightcurve. 

Second, instead of using the polynomial to fit the flux 
versus arclength of the positional drift of the spacecraft, 
we used convolution with a Gaussian kernel, which pro¬ 
vided an improved fit of the correlation compared to the 
polynomial. For each data point the Gaussian convolu¬ 
tion calculates the peak of the Gaussian flux distribution 
of a specified number of data points nearest in arclength 
to a particular data point. At this stage the deviat¬ 
ing data points are typically of astrophysical origin and 
should be rejected to avoid any bias in the SFF calibra¬ 
tion. Trial and error indicated that the flux-versus-drift 
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Figure 2. K2 Campaign 1 short-cadence lightcurve of WASP-85 
before (shown in black) and after the SFF correction (red). Note 
the presence of transits and starspot rotational modulation. The 
corrected lightcurve is shown with an offset of —20 000e"s _1 for 
clarity. 

arclength correlation was best fitted when using a Gaus¬ 
sian kernel with a width of 50 data points and a rejec¬ 
tion threshold of 4 a (see Figure lb). After the Gaussian 
peaks have been calculated for all the data points the 
flux of every data point was replaced with a correspond¬ 
ing Gaussian peak. Figures lc and Id show an example 
lightcurve of a 5-d calibration time interval, before and 
after SFF, respectively. The time interval of 5 d was cho¬ 
sen as the best compromise between including a larger 
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statistical sample of spacecraft drifts and a long-term 
variable drift behaviour. 

With these two modifications we achieved nearly twice 
the precision improvement compared to the original 
kepsff PyKE tool for the short-cadence dataset of 
WASP-85. The median 1-min photometric precision for 
WASP-85 ( K p = 10.2) was 778 parts per million (ppm) 
before SFF and 260 ppm after SFF, an improvement by 
a factor 3. For a comparison, the 1-min photometric pre¬ 
cision of the origin al Ke pler mission for a K p = 10.5 star 
was 202ppm (Still) 2013). 

After removal of all the quality-flagged data points, 
i.e. images affected by the cosmic rays, coarse point¬ 
ing, spacecraft attitude tweak, thruster firing and reac¬ 
tion wheel de-saturation events, and a downlink data gap 
near the middle of Campaign 1, the total number of use¬ 
ful data points was reduced from 120 660 to 108 872. Fig¬ 
ure 2 shows the entire short-cadence lightcurve of WASP- 
85 before and after applying the SFF correction. A total 
of 30 transits are visible. Also visible is the rotational 
modulation caused by the presence of starspots. 


3. ANALYSIS AND RESULTS 
3.1. System parameters 

We obtained the system parameters using the Markov 
Chain Monte Carlo ( MCMC) code presented in [Collier 


Cameron et al. ( 2007|) and further desc ribed in Pollacco 


~t al.]lj2008r ana jAnaerson et ah] ( |2015| ) . In addition to 


the high-quality Iv2 photometry presented herein we used 
the HARPS radial-ve locity measurements of WASP-85A 


from Brown et al. (2014). These measurements are 


not thought to have been significantly contaminated by 
WASP-85B. We corrected the K2 lightcurve for dilution 


by WASP-85B using a flux ratio of 0.50 (see Brown et al. 
( 2014| )). Because the presence of starspot occultations 


may lead to inaccurate d etermination of system parame¬ 
ters (Oshagh et al. 2013) we removed the starspot occul- 
tation events from the measured lightcurve prior to the 
MCMC analysis. Limb darkening was accounted for us¬ 
ing a four-parameter law, with coefficients calculated for 
K epler bandpass and interpolated from the tabulations 
oflSing (2010). 

To refine the ephemeris (but not the other parameters), 
we produced a separate MCMC analysis that also in¬ 
clude d the g round-based photometry presented in Brown 
et al. (|2014 | ), with four-parame ter limb-darkening coeffi- 
cients from [Claret | (2000 2004), as appropriate for differ¬ 


ent bandpasses. This extended the photometric baseline 
from 80 d to 6.5 yr and reduced the uncertainty on the 
period by a factor of 4. 

We present the system parameters in Table 1 and su¬ 
perimpose the corresponding transit model on the K2 
lightcurve in Figure 3. 

There is a discrepancy between the stellar effective 
temperature d erived fro m th e HA RPS spectra (T e ff = 
5685 zb 65 K; Brown et al. 2014) and the value ob¬ 
tained from fitting the transits m our MCMC analy¬ 
sis (T e ff = 6112 =b 27K). The same discrepancy, albeit 
smaller, is present in the discovery paper, which used 
only the ground-b ased tran sits (T e ff = 5910 =b 54 K from 
MCMC; 
etry (eit 


Brown et al. 2014). When fitting the photom- 
rer ground-based or K2) with a transit model 
generated using the spectroscopic T e ff, large systematics 
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Figure 3. Best-fit MCMC model and its residuals from the phase- 
folded K2 lightcurve near the transit. The lightcurve has been 
corrected for dilution by WASP-85B, using a flux ratio of 0.50. 
The lower residual panel corresponds to the effective temperature 
of 6112 K as fitted by our MCMC analysis. The upper residual 
panel serves as a comparison and corresponds to the spectroscopic 
effective temperature of 5685 K which results in large residuals near 
ingress and egress. 

are evident during ingress and egress (see upper residual 
panel in Figure 3), suggestive of a poor modelling of the 
limb darkening. The T e ff discrepancy could be explained 
in part if light from WASP-85B had contaminated the 
spectra of WASP-85A, leading to a lower apparent T e ^ 
for the primary, though any contamination is expected 
to be at a low level. 


3.2. TTV and TDV 

Transit-timing variations (TTVs) and transit-duration 
variations (TDVs) can arise from dynamical interactions 
with additional planets in the system (e.g. Becker et al. 
2015 and citations therein). Typical reported TTV arm 
plitudes for perturbed transiting exoplanets are of the 
order of a few tens of m inutes with periods of a few hun¬ 
dred days (Mazeh et al. 2013]). TDVs on the other hand 
are expected to be observed with smaller amplitudes and 
in phase with TTV. The first unequivocal exoplanet TDV 
detection was measured for system Kepler-88 and indi¬ 
cated a TDV amplitu de of 5 min and a period of 630 d 
(Nesvor ny et al.| |2013). 

We measured the transit-timings and transit-durations 
for WASP-85Ab for each of the transits using the MCMC 
analysis similarly as described in Section 3.1, and show 
the deviations in Figure 4. The detected starspot oc¬ 
cultations have again been removed from the lightcurve 
prior to the analysis as in Section 3.1. The error bars 
in Figure 4 do not include any contribution due to the 
remaining unremoved stellar activity features, which will 
be contributing to the observed scatter of the TTV and 
TDV data points. One data point has been excluded 
from TTV and TDV analysis due to the thruster firing 
event residual that took place at the beginning of ingress 
of the fifth transit near BJD 2456821. 

The hypothesis of the TDV values being zero gives a 
X 2 of 31.8 for 28 degrees of freedom, which means we 
cannot reject it. The y 2 for TTV is 90.6 for 28 degrees 
of freedom which rejects the white noise distribution hy- 
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Table 1 

System parameters from MCMC for WASP-85Ab and its host star 


Parameter 

Symbol 

Value 

Unit 

Transit epoch a 

to 

2456847.472856 ± 0.000014 

BJD 

Orbital period® 

P 

2.65567770 ± 0.00000044 

d 

Area ratio 

(Rp/R*) 2 

0.01870 ± 0.00002 


Transit width 

t T 

n 1 n&i ^+0-00006 
o.iOoib_ 0 00002 

d 

Ingress and egress duration 

t 12, ^34 

fl m Qn^77T0.000008 
U.UIOUO/ ( -0.000016 

d 

Impact parameter 

b 

0 047 +0 ' 005 
u.un/_o 018 


Orbital inclination 

i 

89.691°;“ 

0 

Orbital eccentricity 

e 

0 (adopted) 


Orbital separation 

a 

0.039 ± 0.001 

AU 

Stellar effective temperature^ 5 

T e ff 

6112 ± 27 

K 

Stellar mass 

M* 

1.09 ± 0.08 

M 0 

Stellar radius 

R* 

0.935 ± 0.023 

R© 

Stellar density 

P* 

1.330 ± 0.007 

P© 

Planet equilibrium temperature 

T P 

1452 ± 6 

K 

Planet mass 

M p 

1.265 ± 0.065 

Mjup 

Planet radius 

R p 

1.24 ± 0.03 

Rjup 

Planet density 

Pp 

0.66 ± 0.02 

PJup 


a Epoch and period derived by fitting the photometric datasets from the K2 and all 
the available ground-based observations. 

b This value, obtained from a fit to the transit lightcurve, is discrepant with the 
spectroscopic value (T e ff = 5685 ± 65 K). 
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Figure 4. TTV (upper panel) and TDV (lower panel) for WASP- 
85Ab. The measurements indicate a sinusoidal TTV variability 
with a small semi-amplitude of 7.8 s and a period 14.7 d, likely a sig¬ 
nature of stellar activity rather than actual TTVs. The TDV mea¬ 
surements are consistent with the assumption of constant transit- 
duration. Note that the data point with the largest error bars cor¬ 
responds to the transit during which a reaction wheel de-saturation 
took place near the ingress (see transit 17 in Figure 5). 


pothesis and suggests a sinusoidal TTV signal with a 
semi-amplitude of 7.8 s and a period of 14.7 d. This TTV 
signal, if real, would to the best of our knowledge have 
the smallest detected TTV amplitude to date. On the 
other hand, it is known that stellar activity can induce 


TTV signa ls with semi-amplitudes of up to 200 s (Oshagh 
et al.|2013 ) and that the imperfect removal of some of the 


starspot occultation events, and the effect of the stellar 
rotational modulation, could have resulted in the TTV 
variability. The TTV period of 14.7 d is also compatible 
with the measured rotational period of 15.1 zb 0.6 d (see 
Section 3.4) and it is therefore more likely that the small 
measured TTV signal results from stellar activity. 

The very small TTV amplitude and the absence of sta¬ 


tistically significant TDV indicates that a second, non¬ 
transiting, massive planet is unlikely, unless at a much 
longer period. 

3.3. Starspot detection 

The short 1-min cadence of K2 observations and the re¬ 
stored photometric precision using the SFF enabled us to 
detect individual starspot occult at ions. We show in Fig¬ 
ure 5 the transit lightcurves after subtracting the transit 
model, thus showing only the residuals. We first looked 
for starspots by calculating the y 2 value compared to a 
straight line. This was done by setting a box width equal 
to a phase width of 0.09, and sliding that box in time. 
The box width was determined by examining a few well 
defined starspot occultation events in the system. To 
estimate the threshold of significance we used the same 
box on out-of-transit data. Using this approach we were 
able to detect the 12 most prominent starspot occulta¬ 
tion events above a y 2 ed threshold of 4 (marked with 
asterisks in Figure 5). Below this threshold the residual 
lightcurve proved to be too unstable for this occultation 
detection technique. 

To identify additional occultation events we asked 
four colleagues to mark definite and possible occultation 
events in the residual lightcurves near transits. The or¬ 
der of lightcurves was shuffled to avoid any bias towards 
spots that repeat in consecutive transits. All the events 
which were identified by at least two people as a possi¬ 
ble starspot occultation event are listed in Table 2 and 
marked with red ellipses in Figure 5. We identified a to¬ 
tal of 23 starspot occultations. The measured y 2 ed of the 
individual occultation events are given in Table 2 and 
starspots are marked with different color intensities in 
Figure 5 as a rough indicator of their detectability. 

Blue arrows are added in Figure 5 to mark the po¬ 
sitions of thruster firing events which occasionally intro¬ 
duce artefacts in the lightcurve. Thruster firing artefacts 
are especially present in the first quarter of the observing 
campaign due to the somewhat inconsistent spacecraft 
drifts which might lead to a confusion with occultation 
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Figure 5. Starspot occultations in the model-subtracted lightcurve of WASP-85Ab. Vertical dashed lines specify the extent of the 
transit. Red ellipses mark the 23 identified starspot occultation events. The color intensity denotes the detectability of events and asterisks 
correspond to events detected by our y 2 occultation detection technique. Potential starspot occultation pairs are marked with red lines. 
Blue arrows mark the positions of thruster firing events. Note that transit number 16 lies in the K2 data gap. 


events. Nevertheless, a typical potential thruster firing 
artefact only causes a small step-like deviation in the 
lightcurve and in the out-of-transit WASP-85 lightcurve 
never resembles a starspot-like bump. Therefore, it 
seems unlikely that any of the in-transit starspot-like 
event would be a thruster firing artefact, despite occa¬ 
sional overlaps. 

3.4. Are the starspots recurring? 

The presence of 23 starspots occultation events among 
30 transits indicates that WASP-85A is magnetically ac¬ 
tive, that starspots occur uniformly along the transit 
chord, and that some starspot occultations are poten¬ 
tially occurring in pairs. 

If the stellar rotational period is large compared to 
the orbital period of the planet and if the rotational and 
orbital axes are aligned, it is possible to see an occulta¬ 
tion of the same starspot in consecutive transits. Fig¬ 
ure 5 reveals eight such potential starspot occultation 
pairs at consistent phase shifts. For example, an occul¬ 
tation event at phase 0.013 during transit 11 was likely 
caused by the same starspot that was responsible for the 
occultation event at phase —0.006 in transit 10. All eight 
potential occultation pairs are marked with red lines in 
Figure 5 and their phase shifts and differential stellar 
longitudes are given in Table 2. 


A detection of eight potential occultation pairs (16 
spots) among a total of 23 starspot occultation events 
coupled with the small scatter among measured differ¬ 
ential stellar longitudes indicates that occultation pairs 
were probably caused by occultations of the recurring 
starspots. It is therefore reasonable to suggest that 
planet’s orbital and stellar rotational axis are aligned. 
Using the system parameter values from Section 3.1 we 
estimate that the obliquity angle has to be <14°. From 
the mean differential stellar longitude of the starspot 
pairs of 63.2 d= 2.3° (see Table 2) and known orbital pe¬ 
riod we derive a stellar rotational period of 15.1 =b 0.6 d. 
This is in agreement with the rotational modulation pe¬ 
riod (see Section 3.5) and with the spectroscopically de- 
termin ed rota tional period of 14.2 ±4d for WASP-85A 
( Brown et al. |[2Q14 ). 

3.5. Rotational modulation 

The K2 lightcurve shows a clear rotational modula¬ 
tion (see Figure 2). To estimate the period we used an 
autocorrelation function (ACF). For this task, we first 
removed the transits from the lightcurve and then ap¬ 
plied lOcr clipping. The rotational period is manifested 
as the highest ACF peak near 13.6 d and its multiples 
(see Figure 6). Peaks corresponding to half of this pe¬ 
riod are caused by starspot regions on opposite sides of 
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Table 2 

Starspot occultation events 


Transit 

number 

Med 

Phase 

Stellar 

longitude a [°] 

Pair 

Differential 
longitude [°] 

5 

1.57 

0.0075 

24.9 



6 

7.26 

-0.0077 

-25.6 

6+7 

60.7 ± 3.4 

7 

2.46 

-0.0005 

-1.7 



7 

4.91 

0.0103 

35.1 



10 

6.65 

-0.0059 

-19.4 

10+11 

63.7 ± 4.0 

11 

8.73 

0.0125 

44.3 



13 

5.25 

-0.0037 

-11.9 

13+14 

62.8 ± 4.2 

14 

1.65 

0.0139 

50.9 



19 

5.79 

0.0051 

16.7 



21 

6.21 

0.0004 

1.2 



23 

3.16 

-0.0028 

-9.0 



24 

8.13 

-0.0100 

-33.9 

24+25 

63.8 ± 3.8 

24 

3.06 

0.0062 

20.3 



25 

5.86 

-0.0023 

-7.5 

25+26 

63.3 ± 4.5 

25 

4.21 

0.0089 

29.9 



26 

1.82 

0.0148 

55.8 



27 

2.90 

-0.0136 

-49.4 

27+28 

68.2 ± 4.2 

28 

3.10 

0.0058 

18.8 



29 

5.63 

-0.0028 

-9.0 

29+30 

62.1 ± 4.3 

29 

2.31 

0.0031 

10.0 



30 

4.95 

-0.0064 

-20.9 

30+31 

61.3 ± 3.8 

30 

2.29 

0.0143 

53.1 



31 

1.81 

0.0116 

40.4 




a Stellar longitude of zero corresponds to the meridian that runs 
through the center of the stellar disc. 


the star. The 80-d time-span of the K2 observations cov¬ 
ers five cycles of the 13.6-d periodicity. Following; the 
procedure demonstrated by McQuillan et al. (2013) we 
determined the rotational period as a median of the in¬ 
tervals At between the ACF peaks associated with the 
rotational period. We conservatively estimated the error 
as the FWHM of the highest ACF peak. 

The resulting rotational period of 13.6 ±1.6 d is in 
agreement with the spectroscopically determined rota¬ 
tional period for WASP-85A, 14.2 ± 4 d, if the star’s ro¬ 
tational axis is orthogonal to the line of sight. But, the 
value is also in agreement with the estimated rotational 
period for WASP-85B, 11.7±4d (Brown et al. 2014|). 
Therefore, it is uncertain which of the two binary com¬ 
ponents is causing the observed rotational modulation. 
However, because WASP-85A is active (see Figure 5), 
it is reasonable to assume that the rotational modula¬ 
tion is caused by the rotation of the spotted surface 
of WASP-85A. Because rotational modulations depend 
strongly on starspots coming and going, which causes 
the phase shifts, the period derived from it is therefore 
less accurate and less reliable than the period determined 
via starspot occultations (see Section 3.4). 

We used S pot O scillat ion And Planet 2.0 (SOAP 2.0, 


Dumusque et al. (2014)) to fit the starspot rotational 
modulation of WASP-85 and found that the three char¬ 
acteristic “M”-shaped modulation features can be well 
fitted with two main starspot regions on the opposite 
sides of the rotating stellar surface. However, the change 
in the lightcurve over time shows that the starspot re¬ 
gions are changing over the K2 observations. 


3.6. Phase-curve variations 

We folded the K2 lightcurve on the orbital period of the 
planet. In order to do this we first removed the rotational 



Figure 6. Autocorrelation function of the K2 Campaign 1 short- 
cadence lightcurve of WASP-85. Marked with P is the first and 
highest peak corresponding to t3.6d rotational period. Additional 
four peaks are visible at multiples of this period. Peaks corre¬ 
sponding to half period are caused by another starspot region on 
the opposite side of the star. 


modulation using our two-region model of the starspots. 
However, since this model is only an approximation, the 
folded lightcurve is still dominated by the residuals of the 
starspot modulation, which do not average out owing to 
the limited length of the K2 observations (Figure 7). 

We did not detect the secondary eclipse nor any other 
phase-curve variations caused by the orbiting exoplanet. 
Using the system parameters from Table 1 the ex¬ 
pected semi-amplitudes of ellipsoidal, Doppler, and re¬ 
flection variations are 3.7, 2 .1, an d (for an albedo of 
0.1) 15.4 ppm, respectively (Pfahl et al.|2008|). We con¬ 
servatively estimate an upper limit - ©! 50 ppm for the 
semi-amplitude of variations over the planet’s phase- 
curve. Figure 7 shows a comparison between the mea¬ 
sured phase-curve of WASP-85 and the simulated reflec¬ 
tion modulation with a semi-amplitude of 50 ppm. The 
inability to detect the reflection modulation and the sec¬ 
ondary eclipse in the measured K2 phase-curve of WASP- 
85 indicates that the planet’s visual geometric albedo is 
lower than ^0.3. 


4. CONCLUSION AND DISCUSSION 

We have modified the PyKE self-flat-fielding method to 
nearly double its performance in restoring the degraded 
short-cadence photometric precision due to the pointing 
drift of the K2 spacecraft. We find similar syste m param¬ 
eters to the discovery paper Brown et al. (2014), measure 
a rotational modulation period of 13.6 d, and reconstruct 
the “M”-shaped photometric modulations with two dy¬ 
namic starspot regions on either side of a magnetically 
active star. The presence of starspot occultations indi¬ 
cates that the host star WASP-85A is active and sug¬ 
gests that the stellar rotation and the planet’s orbital 
axis are aligned. The phase shifts of potential occul¬ 
tation pairs correspond to a stellar rotational period of 
15.1 ± 0.6 d. The agreement between the stellar rota¬ 
tion and transit-timing variation period suggests that its 
small semi-amplitude of 8 s is a consequence of stellar ac¬ 
tivity rather than an inter-planet gravitational interac¬ 
tion. Measured transit-duration variability is consistent 
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Figure 7. Measured phase-curve of WASP-85, binned to 200 bins, 
(shown in black) and a simulated reflection phase-curve modulation 
with a semi-amplitude of 50ppm (red). The inability to see the 
reflection modulation and secondary eclipse in the measured phase- 
curve of similar amplitude as shown in the simulated phase-curve 
suggests that the actual reflection semi-amplitude is lower than 
50 ppm. 


with white noise distribution around zero. We do not 
detect any phase-curve modulations nor the secondary 
eclipse of the planet. 

Starspots provide one method of constraining the sky- 
projected obliquity of the host star, while the Rossiter- 
McLaughin (R-M) effect provides an independent esti¬ 
mate. 

For example, in WASP-4 the R-M effect i mplies a pro- 
jected alignment angle of A = — 4^34 (Triaud et al. 


2010), while Sanchis-Ojeda et al. (2011) report that 
bn two occasions the same starspot was occulted, from 
which they constrain the alignment to A = — l 0 !^- Sim¬ 
ilarly in WASP-6 the R-M angle is A = (Gillon 


et al.|2009 ), and Tregloan-Reed et al. (2015) suggest that 


the same starspot was occulted in two transits and re¬ 
port A = 7.2 ± 3.7°. In WASP -19 the R-M angle is 
A = 4 . 6 ± 5. 2° ( Hellier et al.||2011 ), while Tregloan-Reed 


et al. (2013) saw the same starspot occulted on consecu- 
tive nights and report A = 1.0 db 1.2°. 

HAT-P -11 is misaligned, with an R-M angle of A = 
103°^io (Winn et al. 2010j). The mis alignment is con- 
firmed by th e analysis of starspots by Sanchis -Ojeda fe 


Winn (|2011 ). Similarly, for Kepler-63 Sanchis-Ojeda 


et al. 


(2013) report A = — llO 0 ^ from analysing both 
starspot”and the R-M effect. 

There is currently no measurement of the R-M effect 
in our system, WASP-85A, but we predict that it would 
confirm that the system is aligned. Qatar-2 is a similar 


case, where starspots indicate alignment (Mancini et al 


2014) though this is not yet confirmed by an R-M mea¬ 


surement. 

There are no reported cases where the R-M angle 
conflicts with the analysis of starspots, and thus the 
starspots are a reliable way of measuring alignment in 
systems too faint for R-M measurements. Examples 
are Kepler-30 ( V = 15.6) where occultations of starspots 
by both Kepler-30c a nd Kepler-30d indicate alignment 
(Sanchis-Ojeda et al. 2012|, and Kepler-1 7 (V = 14.3) 


where starspots also indicate alignment (Desert et al. 


2011 ). 

The benefit of ongoing K2 observations of known 
WASP planets is that the systems are generally bright 


enough to enable good measurements of the alignment 
angle, A, the stellar rotational period, the stellar radius, 
and the star’s projected rotational broadening, vsm i. 
Given all of these we obtain, not only the projected align¬ 


ment A, but the actual alignment angle ip (see, e.g., Winn 
et al 2005). 


Adopting the rotational period of 15.1 ± 0.6 d, and 
combining with a stellar radiu s of 0 . 94 R^ (Ta ble 1) and 
a v sin i of 3.41 d= 0.89 kms -1 (Brown et al. 2014) we find 
that the star’s rotational axis must be at >50 u to the line 
of sight. This leads to the constraint that 'ip < 41°. 

The measured effective temperature and spin-orbit 
alignment of WASP-85 system are consistent with the 
current sample of predominantly aligned systems with 
stars below the characteristic effective temperature of 
6250 K. 
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